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Evolution of minimum-bias parton fragmentation in nuclear collisions 
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Minimum-bias fragment distributions (FDs) are calculated by folding a power-law parton energy 
spectrum with parametrized fragmentation functions (FFs) derived from e^-e~ and p-p collisions. 
Changes in FFs due to parton "energy loss" or "medium modification" are modeled by ahering FF 
parametrizations consistent with rescaling QCD splitting functions. The common parton spectrum 
is constrained by comparison with a p-p pt spectrum hard component. In- vacuum and in-medium 
FDs are compared with spectrum hard components from 200 GeV Au-Au coUisions for several 
centralities. The reference for aU nuclear coUisions is the FD derived from in-vacuum e'^-e~ FFs. 
The hard component for p-p and peripheral Au-Au collisions is found to be strongly suppressed for 
smaller fragment momenta, consistent with the FD derived from in-vacuum p-p FFs. At a particular 
centrality the Au-Au hard component transitions to enhancement at smaller momenta and suppres- 
sion at larger momenta, consistent with FDs derived from in-medium e'^-e~ FFs. Fragmentation 
systematics suggest that QCD color connections change dramatically in more-central A-A collisions. 
Spectrum systematics are inconsistent with saturation-scale arguments and parton thermalization. 

PACS numbers; 12.38.Qk, 13.87.Fh, 25.75.Ag, 25.75.Bli, 25.75. Ld, 25.75.Nq 
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I. INTRODUCTION 



RHIC collisions are conventionally described in terms 
of two major themes: hydrodynamic evolution of a 
thermalized bulk medium [l|, |3, H, 01 and energy 
loss of energetic^ partons in that medium via gluon 
bremsstrahlung [3]. Medium dynamics and properties 
and parton specific energy loss relating to "tomography" 
of the medium are the principal analysis goals [6|, |7|. 
Analysis methods tend to favor those goals: Methods 
directed toward a bulk medium tend to suppress low-pt 
features of parton fragmentation, and methods applied 
to high-pt jet analysis also tend to suppress structure at 
smaller pt @] . 

Recent physical-model-independent studies of spec- 
trum and correlation structure have revealed interesting 
new aspects of RHIC collisions. Analysis of number and 
Pt angular correlations led to unanticipated structure in 
the final state, subsequently identified with parton frag- 
mentation in the form of minijets [lol . [ill , [l^ . Il3| |. 
Two-component analysis of p-p and Au-Au spectra re- 
vealed a corresponding hard component, a minimum-bias 
fragment distribution associated with minijets which sug- 
gested that jet phenomena extend down to 0.1 GeV/c \§, 

In this analysis new aspects of spectra and correla- 
tions in p-p and Au-Au collisions arc combined with com- 
plete representations of fragmentation functions from e^- 
e~ and p-p [3, [13 collisions to reveal the systematic 
evolution of parton fragmentation with Au-Au centrality. 
Accurately parametrized FFs are combined with a power- 
law parton spectrum to produce fragment distributions 
which can be compared quantitatively with hard compo- 
nents derived from pt spectra in nuclear collisions. The 
observed FD evolution reveals surprising new features of 
parton fragmentation in p-p and A-A collisions. 



II. MINIJETS 

Minijets dominate the transverse dynamics of nuclear 
collisions above ^/SNN ~ 15 GeV. They have an ex- 
perimental and theoretical history of more than twenty 
years. The term "minijets" can be applied collectively to 
all hadron fragments from the minimum-bias scattered- 
parton spectrum averaged over a given A-A or N-N event 
ensemble. Because the parton pt spectrum is rapidly 
varying (~ 1/Pt)i the minimum-bias spectrum is nearly 
monoenergetic, peaked at an effective termination or cut- 
off energy near 3 GeV [H, [l9, 20]. The term "minijets" 
then refers experimentally to jets localized near the cutoff 
energy. 

Minijets manifest as both minimum-bias jet corre- 
lations [2l[ and as the hard component of the two- 
component spectrum model [T^ . They provide unbi- 
ased access to fragment distribution structure down to a 
small cutoff energy for scattered partons (those fragment- 
ing to charged hadrons) and to the smallest detectable 
fragment momenta 0.1 GeV/c). Because they are 
large enough to observe accurately but small enough to 
respond fully to any QCD medium, minijets serve as 
Brownian probes of QCD in nuclear collisions [l^l . 

The minijet concept emerged experimentally at the 
SPPS from a UAl analysis of Et structure down to 
small integrated Et [l^. The analysis determined that 
Et "clusters" (minijets) are distributed according to the 
expected pQCD power-law parton spectrum down to 5 
GeV. Azimuth correlations between clusters exhibited a 
peak at tt radians expected for back-to-back scattering 
of initial-state partons. The 5 GeV Et cutoff was later 
related to a 3-4 GeV parton energy equivalent [20| . 

Corresponding minijet structure was observed in two- 
particle correlations from 200 GeV p-p collisions [2]| . 
Angular correlations with no "jet" pt conditions exhibit 
just the structure expected from pQCD jets: a narrow in- 
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trajet same-side peak at the angular origin (parton frag- 
mentation) with most-probable pt ~ 1 GeV/c and an 
interjet away-side ridge at tt radians (back-to-back par- 
ton scattering). As noted by UAl, there is no dividing 
line between conventional high-pt "jets" and "minijets." 

Several theoretical treatments have identified experi- 
mental minijets with parton scattering and fragmenta- 
tion. Minijet production was calculated perturbatively 
for anticipated RHIC U-U collisions based on the UAl 
minijets: "The observed [minijet] rate is in agreement 
with [p]QCD and is quite large" "Semihard par- 

ton interactions [as in the UAl minijet analysis] appear 
to play an important role in high-energy hadronic scat- 
tering [y^ > 10 GeV]" A "...theoretical cutoff 
of p™*" 3 GeV seems to describe the observed total 
minijet cross section with Ei^* {E^°'"') > 5 GeV" and pro- 
duces a minijet total cross section in agreement with UAl 
data ^3]- "There is an increasing amount of evidence 
that the perturbative domain of QCD extends down to 
[parton] momenta of the order of 1 GeV/c" [HjUHl- 

Minijet thermalization is of central interest as the ba- 
sis for QGP formation and hydrodynamic flows. "Mini- 
jets... will be reprocessed by the system and not emerge 
from it" [l^ . In [2^ the thermalization time is estimated 
as 4-5 fm/c with T ~ 200 GeV. One basis for claims of 
thermalization is the assumption that partons (gluons) 
propagate in a gas of bare gluons, which can be strongly 
questioned given recent minijet-related results at RHIC 
(e.g., the present analysis and [ol. [lol. HI. [isl. [2]|). 

The HIJING Monte Carlo was developed specifically 
to study the role of minijets in p-p and A-A collisions: 
"We emphasize the effects due to multiple mini-jet pro- 
duction at collider energies" [27,, 28] . HIJING p-p corre- 
lations quantitatively match minijet correlations (same- 
side amplitude, widths, away-side ridge) measured in p-p 
collisions |21]. HIJING predictions with "jet quenching" 
disabled are consistent with a Glauber linear superposi- 
tion reference for A-A collisions [10, E3l ■ 



III. ANALYSIS METHOD 

In this analysis minijets manifested as pt-spectrum 
hard components are the object of study. Experimen- 
tal hard components are modeled with fragment distri- 
butions calculated by folding parton spectra with var- 
ious fragmentation-function ensembles. From the com- 
parisons parton spectrum parameters and modifications 
to fragmentation functions in more-central Au-Au colli- 
sions are inferred. 

A study of the charge-multiplicity Uch dependence of 
p-p pt spectra revealed two components with fixed func- 
tional forms independent of n^h, denoted soft and hard 
components and later interpreted in terms of longitudinal 
projectile fragmentation (soft) and transverse scattered- 
parton fragmentation (hard) [T^ . Separation into two 
components was based on a Taylor expansion of spec- 
trum structure on event Uch with no physical model im- 



posed. The hard-component fragment distribution or FD 
was identified with minijet angular correlations, having 
all the characteristics of jet correlations but with nojet- 
specific cuts imposed (i.e., minimum-bias jets) |2l| . 

Analysis of e'^-e~ (e-e) fragmentation functions (FFs) 
from LEP and HERA led to a complete characteriza- 
tion of FFs for all parton energy scales in terms of beta 
distributions on normalized rapidity u. FFs were rep- 
resented to their statistical limits for all fragment mo- 
menta (in contrast to the limitations of conventional 
pQCD parametrizations such as the MLLA) [l^. The 
most important achievement was accurate representation 
of FFs near the smallest fragment and parton momenta. 
Comparison with p-p (p-p) FFs indicated substantial sys- 
tematic e-e vs p-p differences for smaller fragment mo- 
menta not revealed by conventional data plots on mo- 
mentum fraction Xp — p/pjet or its logarithmic equiva- 
lent = ln(l/xp). 

In this analysis e-e and p-p FFs are folded with a par- 
ton spectrum model to produce fragment distributions 
to be compared with measured spectrum hard compo- 
nents. For example, a measured FF ensemble from p-p 
collisions is folded with the spectrum model to produce 
a calculated FD. Comparison of the FD with a measured 
p-p spectrum hard component determines the spectrum- 
model cutoff energy and QCD power-law exponent. The 
parton spectrum agrees quantitatively with pQCD pre- 
dictions and is comparable with a UAl measurement of 
the differential jet cross section based on Et clusters. 

The question then arises which FFs should be used to 
calculate FDs for various nuclear collision conditions. In 
the initial part of the analysis it is assumed that the un- 
derlying parton spectrum remains unchanged for a given 
final-state hadron acceptance (e.g., pt > 0.15 GeV/c, 
\ri\ < 1, 27r azimuth). Manifestations of "jet quenching" 
or parton "energy loss" are modeled via FF modifications 
accessible for the first time over the full pt acceptance. 
A scheme for FF medium modification in '29| is found to 
be particularly relevant to data. 

Spectrum ratio measures and direct comparisons of 
hard components with calculated FDs reveal the rela- 
tion between systematic FF modifications and Au-Au 
centrality. Ironically, the most significant fragmentation 
changes occur below ~ 2 GeV/c where most fragments 
appear but where comparisons with pQCD are typically 
deferred in favor of hydrodynamic descriptions. In partic- 
ular, a direct correspondence has emerged between dra- 
matic fragmentation modifications noted in this analysis 
and the recently-observed sharp transition in minijet cor- 
relation systematics at a specific Au-Au centrality [l3] • 



IV. P-P TWO-COMPONENT SPECTRA 



The two-component model of p-p spectra [14| was the 
starting point for the differential fragmentation analy- 
sis described in this paper. The two-component model 
emerged from a Taylor-series expansion of spectra on 
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observed charge multiplicity rich in one unit of rj (~ 
dhch/drj) and was not motivated by a particular physical 
model. The spectrum components (Taylor series coef- 
ficients) were subsequently interpreted physically in the 
context of correlation analysis and by analogy with par- 
ton fragmentation systematics at larger energy scales. 



A. Two-component spectrum model 

The two-component model applies to two-particle cor- 
relations and to their ID projections, the pt or yt spectra. 
The two-component spectrum model for p-p collisions 
sorted according to event-multiplicity index rich is 
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where n,j; is integrated over one unit of rj (i.e., nx/2'K ^ 
d'^nx/dr]d(/)), soft component S'o(yt) is the Taylor series 
"constant," and hard component -ffo(2/t) is the coefficient 
of the term linear in rich, both normalized to unit integral. 
For comparisons with A-A spectrum data we define Spp = 
(l/yt) dris/dyt with reference model ris So and similarly 
for Hpp <-> rihHo. The two-term Taylor series exhausts 
all significant p-p spectrum structure. 

Fig. [1] (left panel) shows spectra for ten multiplicity 
classes from 200 GeV non-single diffractive (NSD) p-p 
collisions The asymptotic limit for hch (dash- 
dotted curve) is 5*0. The spectra are normalized by the 
soft-component multiplicity = rich/i^ + cxfi'ch), where 
a ~ 0.01 and hch ~ 'n-ch/'i is the observed rich resulting 
from incomplete pt acceptance and tracking inefficiencies. 




FIG. 1: Left: Spectra for ten multiplicity classes [1,11.5] of 
200 GeV NSD p-p collisions The dash-dotted curve is 

the spectrum soft component So{yt) defined as the limiting 
case for hch — > 0. Right: The two-component (soft-|-hard) 
model of p-p spectra. Hard component Ho{yt) is a Gaussian 
with QCD power-law tail 

Fig. [T] (right panel) shows the two-component algebraic 
model Eq. ([T]) with unit-normal model functions 5*0 and 
Ho defined in ^ [ij]. The hard-component spectrum 
contribution rih/ris scales as ahch- Factor a = 0.01 is 
the average value for most hch classes. The factor drops 
to 0.0055 for hch = 1- The spectrum data in the left 
panel are described to the statistical limits. 



B. p-p spectrum hard component 



Figure [5] (left panel) shows p-p hard components in 
the form Hpp/ris for ten multiplicity classes obtained by 
subtracting fixed soft component So from the ten NSD 
p-p spectra normalized by soft multiplicity rig . The hard- 
component shape is independent of multiplicity and de- 
scribed approximately by a Gaussian. The amplitude is 
approximately proportional to hch- From the figure Ho 
coefficient rih/ris — ahch is inferred, with a ^ 0.01 [l3|. 




FIG. 2: Left: Spectrum hard components Hpp{yt,hch) (solid 
points) for ten multiplicity classes of 200 GeV p-p colli- 
sions [Ij]. A common Gaussian model function describes the 
data well except for parts of the lowest multiplicities. Right: 
Hpp data from the left panel normalized to NSD p-p collisions 
by factor nh{l -25) /rih (hch) demonstrating the common form. 
A Gaussian (dotted curve) and Gaussian with power-law tail 
(dash-dotted curve) are compared to the data. 



Figure [2] (right panel) shows ten hard components Hpp 
from the left panel scaled by factors nh{l.25)/nh{hch) to 
reveal a common shape representing the mean hard com- 
ponent for NSD p-p collisions. The dash-dotted curve is 
0.02 i7o [0.02 - (a = 0.007) {hch = 1-25) (n^ = 2.5) [H], 
with Ho defined as a Gaussian plus exponential tail on 
transverse rapidity yt (rih Hq is relabeled below as ref- 
erence Hgg)- The exponential tail represents the QCD 
power law oc exp[(nQCD — 2) yt], where the —2 

results from the pt — > yt Jacobian [8|] . The dotted curve, 
a Gaussian with no QCD tail [l3], is inconsistent with 
data at larger yt- The spectrum hard component is in- 
terpreted as a minimum-bias fragment distribution dom- 
inated by "minijets" — jets from those partons (gluons) 
with at least the minimum energy required to produce 
charge-neutral combinations of charged hadrons. 



V. FRAGMENTATION FUNCTIONS 

e-e and p-p FFs for inclusive fragments and inclusive 
partons have been parametrized accurately over the full 
{y,ymax) ^ ixp,Q'^) region relevant to nuclear colli- 
sions. The parametrizations permit comprehensive tests 
of pQCD in relation to nuclear collision data. 
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A. Accurate FF parametrizations 

Fragmentation functions provide direct experimental 
access to the parton-hadron interface of QCD. At RHIC 
energies nuclear collisions are dominated by parton scat- 
tering and fragmentation. For full understanding of col- 
lision dynamics FFs should be described over the entire 
fragment distribution and over all parton energies rele- 
vant to nuclear collisions. 

pQCD theory emphasizes the parton splitting cascade 
and DGLAP evolution of FFs. The 10% most-energetic 
fragments are well-described (e.g., [l^l). FFs are conven- 
tionally represented by semilog plots of D(xp, Q^) on mo- 
mentum fraction Xp = p/pparton, which compress and ob- 
scure details at small Xp (including most of the fragment 
distribution) or on ^p = ln(l/a;p) for which the lower lim- 
its of FFs on p are not well-defined. In contrast, rapidity 
y = ln[(_E +p)/mT^] is well-defined as fragment p ^ 
iy ~^ p/fTT-Tr), which is essential for study of FFs in nu- 
clear collisions. The small-Xp (small-y) region dominates 
nuclear collisions, both by driving collision dynamics and 
by providing diagnostic evidence in the final state. 

Two goals should be distinguished: 1) phenomenologi- 
cal descriptions of FFs which can provide simple and ac- 
curate representations of data over the large kinematic in- 
tervals required for comprehensive study of nuclear colli- 
sions; 2) theoretical descriptions of FFs via pQCD, which 
may be limited by the current state of theory. 



B. e'''-e fragmentation functions 

e+-e^ light-quark (uds) and gluon fragmentation func- 
tions are well-described above energy scale (dijet en- 
ergy) Q ~ 10 GeV by a two-parameter beta distribu- 
tion f3{u]p, q) on scaled rapidity u Q^. Parameters {p, q) 
vary slowly and linearly with Q above 10 GeV and can 
be extrapolated reasonably well down to Q ~ 4 GeV 
based on dijet multiplicity data. Most of the FF "scal- 
ing violations" described by the DGLAP equations re- 
sult from self-similar variations of FFs with energy scale, 
which can be absorbed into dijet multiplicity 2 n{Q) and 
scaled rapidity u = {y - ymin)/{ymax - Vmin), with 
ymax = ln(Q/m^) and fixed j/mm inferred from the sys- 
tematics of measured e~^-e~ FFs. Dijet multiplicities are 
determined by P{u;p, q) according to an energy sum rule. 

Fig. [3] (left panel) shows measured FFs for three energy 
scales from HERA/LEP [Hi HI]. The 2 in the axis label 
indicates that these are dijet rich densities. The vertical 
lines at right denote ymax values. The curves are deter- 
mined by the {p, q) parametrization with ymin ^ 0.35 
{pt ~ 0.05 GeV/c, left vertical line) and describe data 
to their error limits over the entire fragment momentum 
range. 

Fig. [3] (right panel) shows the FF ensemble (inclusive 
light quarks fragment to inclusive hadrons) vs energy 
scale Q as a surface plot. The dashed curve is the lo- 
cus of modes — the maximum points of the FFs. Between 
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FIG. 3: Left: Measured e'^-e~ fragmentation functions (sym- 
bols) for three energy scales (dijet energies) [3ll. |3^. Curves 
through data are from a universal parametrization based on 
the beta distribution [l^. Right: (Color online) Surface 
plot of the universal FF parametrization showing the locus 
of modes (dashed curve). 



the dash-dotted lines the system is determined by fiducial 
FF data (i.e., exceptional quality and range T5]) Between 
the dash-dotted and dotted lines the parametrization is 
constrained only by dijet multiplicities. 

Systematic trends can be extrapolated to the left of 
the left dotted line. At the end of a splitting cascade 
partons evolve to color-singiet hadrons 1-to-l with a sin- 
gle hadron in the parton "fragmentation function" and 
y fragment = yparton (diagonal fine). Approach of the 
locus of modes to the diagonal therefore leads to local 
parton-hadron duality (LPHD) [s^. At lower energy 
scales the QCD density of states is small and hadron 
resonances dominate. The QCD splitting cascade tran- 
sitions to a resonance cascade terminating in detected 
hadrons. 
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FIG. 4: Left: parametrized e"'"-e~ fragmentation functions for 
five dijet energies. Right: Integrated dijet multiplicities for 
in-vacuum FFs in the left panel (solid points) and for FFs 
modified according to a model of "energy loss" or medium 
modification |29| in central Au-Au collisions (open circles). 



Figure [3] (left panel) shows parametrized beta FFs for 
five energy scales. The 6- GeV scale is relevant to the 
minijet spectrum and fragment distributions from this 
analysis. Such curves provide a complete description of 
e-e FFs at energy scales relevant to nuclear collisions. 
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Figure |4] (right panel) shows light-quark dijet multi- 
plicity systematics from the same beta parametrization. 
The solid points correspond to the FFs in the left panel. 
The open circles represent multiplicities from medium 
modification of those FFs in central Au-Au collisions at 
200 GeV, as described in Sec. IVIIIBI The "in-medium" 
shift of FFs to smaller fragment momenta requires more 
fragments to satisfy energy conservation. The system- 
atics of quark and gluon jets coincide for energy scales 
Q = 2Ejet < 8 GeV. Quark-gluon differences at larger 
energy scales are less important for gluon-dominated 
minimum-bias fragmentation in nuclear collisions. 



C. p-p fragmentation functions 

Figure El (le ft panel) shows FF data from p-p collisions 
at FNAL The plotted points are samples from the 
full data set. Jets are integrated within a cone half- angle 
of 0.47 radians. The solid curves guide the eye. There 
is a significant systematic difference between p-p and e- 
e FFs. The dotted line represents the lower limit for 
e-e FFs. The systematic gap for all parton energies is 
apparent — ymin for p-p colhsions is ~ 1.5 (0.3 GeV/c) 
instead of 0.35 (0.05 GeV/c). The curve labeled MB is 
hard-component reference Hqq from p-p collisions [3], 
comparable to the 6 GeV FF curve in Fig. [5] (left panel). 
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FIG. 5: Left: Measured fragmentation functions (samples) 
from p-p collisions at several energies (symbols) [iq ]. Solid 
curves guide the eye. The dotted line represents contrasting 
e^-e~ FF systematics. Right: (Color online) Surface plot of 
the universal e'''-e~ FF beta parametrization modified with a 
common cutoff factor to describe p-p FFs. 

Figure [5] (right panel) shows a surface plot of the p-p 
FF ensemble. For a systematic representation of p-p FFs 
the e-e FF beta parametrization has been modified by 
adding a cutoff factor 



ffcut(y) = tanh{(y -yo)/(,y} y > ya, 



(2) 



with 2/0 ~ ~ 1-6 determined by the CDF FF data [l6|. 
The modified e-e FFs have not been rescaled to recover 
the initial parton energy. The cutoff function represents 
real fragment and energy loss from p-p relative to e-e 
FFs. 



Figure [BKleft panel) shows e-e beta FFs for five parton 
energies [1^ modified by the gcut factor. Also plotted 
are more-recent CDF FF data for dijet energies 101 and 
216 GeV [l3| demonstrating the correspondence. The 
CDF FFs also reveal a systematic amplitude saturation 
or suppression at larger parton energies compared to LEP 
systematics, evident also in Fig. [5] (left panel). In Fig.lH] 
(left panel) the 216 GeV p-p data fall well below the 216 
GeV LEP expectation (sohd curve). The 101 GeV p- 
p data fall below the 101 GeV e-e parametrization to a 
lesser degree, and mainly for smaller fragment rapidities. 
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FIG. 6: Left: Parametrized p-p fragmentation functions for 
five dijet energies (solid curves). The FF data (samples) at 
two energies (symbols) reveal significant suppression rel- 
ative to the e^-e~ beta parametrization. Right: Calculated 
dijet multiplicities (solid points) for in-vacuum p-p FFs (solid 
curves in left panel) showing significant reductions from the 
in-vacuum e^-e~ trend (solid curve). Open triangles [13] and 
open circles [35| are CDF measured p-p dijet multiplicities. 

Figure [S] (right panel) shows multiplicity systematics 
(solid points) for p-p (i.e., modified e-e) FFs. The solid 
curve represents unmodified e-e FFs. There is substantial 
reduction of p-p FF multiplicities due to the cutoff. Also 
plotted are CDF FF multiplicities from reconstructed jets 
(open triangles 34] and open circles [ssj). 

At 100 GeV dijet energy {Ejet = 50 GeV) the FF mul- 
tiplicity in p-p collisions is reduced by ~ 6 relative to 
e-e FFs, with missing-fragment mean pt ^ 0.4 GeV/c. 
The corresponding ~ 2.5 GeV missing energy represents 
a small fraction of the total jet energy (possibly within a 
calorimeter calibration error). But the 30% rich reduction 
could have a major impact on the description of nuclear 
collisions. At smaller energy scales the fractional multi- 
plicity reduction is much larger (e.g., 70% for Ej^t = 3 
GeV). We return to that important issue in Sec. IXI Bl 



VI. FRAGMENT DISTRIBUTIONS 

Whereas a fragmentation function (FF) is conditional 
on a specific parton energy, a fragment distribution (FD) 
is the hadron distribution associated with a minimum- 
bias parton spectrum — the folding of an FF ensemble 
with the parton spectrum. The hard component from p-p 
spectra can be interpreted as an FD ^14] , consistent with 



6 



p-p correlation systematics [2i|. This analysis provides 
further support for the FD interpretation. 

NLO "fragmentation functions" [sHi combine a pQCD 
parton spectrum with a theoretical parametrization of in- 
vacuum e-e FFs, e.g. as measured at LEP/HERA |^3l|,l3^. 
Two questions arise: 1) Is the theory description of FFs 
adequate over the entire fragment momentum range and 
parton energy range relevant to nuclear collisions? 2) Are 
e-e FFs appropriate for p-p collisions — is the assumption 
of FF universality inherent in the QCD factorization the- 
orem relevant? The reply to 1) is currently no. The reply 
to 2) depends on context, as revealed by this analysis. 



represents in this analysis the minimum parton momen- 
tum which leads to detectable charged hadrons as neutral 
pairs (i.e., local charge ordering [SSfV 
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A. The pQCD folding integral 



The folding integral used to obtain FDs in this analysis 



d'^rih ^ ejSr], A77) 
dy dri (Jnsd A77 7o 



dyniax ^xxi^y^ ymax) 



dadi 



dijet 



dy„ 



(3) 



where D:xx{y,ymax) is the dijet FF ensemble from a 
source collision system (xx ~ e-e, p-p, AA, in-medium 
or in- vacuum), and da dijet /dymax is the minimum-bias 
parton spectrum. Spectrum hard component d^nh/dy dr] 
as defined represents the fragment yield from scattered 
parton pairs into one unit of rj. Efficiency factor e ~ 0.5 
(for a single dijet and one unit of rj) includes the prob- 
ability that the second jet also falls within 77 acceptance 
St] and accounts for losses from jets near the acceptance 
boundary. A77 ~ 5 is the effective 47r 77 interval for scat- 
tered partons. Further details are given in Sec. IVIII 



B. Parton spectrum model 

The effective parton spectrum for charged hadron frag- 
ments from p-p collisions can be estimated by folding a 
pQCD power-law parton spectrum hypothesis with trial 
FFs from p-p and e-e collisions and comparing the re- 
sulting FDs with the measured p-p spectrum hard com- 
ponent (interpreted as a fragment distribution). 

A model for the parton pt spectrum resulting from 
minimum-bias scattering into an rj acceptance near pro- 
jectile mid-rapidity can be parametrized as 



1 dad 



ijet 
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Pt dpt 



riQCD ' 
Pt 



(4) 



which defines exponent uqcd- The equivalent jet spec- 
trum on ymax = ln(2pt/m^) is 



da, 



dijet 



dy„ 



P't 



Pt riQcn 
Pt 



= ^Vma^ exp{-(nQCD - 2)y„iaa;}, 

where p^ is the Jacobian factor for pt ymax and j/n 
ycutj the spectrum cutoff. The cutoff factor 



(5) 
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FIG. 7: Dijet (parton-pair) transverse energy spectra on ra- 
pidity ymax ~ \ri{2 Ej St/ m-^) plotted in semilog (left) and lin- 
ear (right) formats. The solid curves are determined by a mea- 
sured p-p spectrum hard component. The dash-dotted curves 
illustrate reduction of the cutoff energy inferred for central 
Au-Au collisions. The bold dotted curve labeled pQCD in 
the left panel is discussed in Sec. Ixra The light dotted ex- 
trapolation down to 1 GeV corresponds to a saturation-scale 
cutoff estimate CSec. IXIII"G|) . 

Fig. [7] (semilog and linear formats) shows a parton 
spectrum inferred from the p-p spectrum hard compo- 
nent in the next subsection, ycut is well-defined by the p- 
p hard component, and uqcd is defined by Au-Au spec- 
trum hard components extending to larger yt . Width pa- 
rameter ^cut affects details of FDs below the maximum 
(mode). Fixed value <^c«t = 0.1 produces FD shapes con- 
sistent with ^cut = but avoids a discontinuity. 

For a given value of jet cross section adijet the coeffi- 
cient is Ay^^^ = adijet {riQCD - 2) exp{(nQCD - 2) ycut}- 
For nominal values uqcd = 7.5 and ycut = 3.75 {Ecut ~ 
3 GeV) adijet ~ 2.5 mb is adjusted to match FD data at 
larger y (power-law tail) , defining a fixed value of Ay^^_^ . 
ycut is then adjusted to fit FD structure near the mode. 
The actual jet cross section varies strongly with ycut as 



dad^ 



[ijet 



^ dijet 



^{nqcD - 2) dycut- 



(7) 



A 0.1 reduction in j/cwt, e.g. from 3.75 to 3.65, (10% rela- 
tive reduction in Ecut) leads to a 55% increase in the jet 
cross section (cf. Fig. [11]- left panel). 



C. Fragment distributions from p-p and e-e FFs 

Specific FF ensembles from LEP/HERA e-e and FNAL 
p-p collisions can be combined with the parametrized 
parton spectrum to produce FDs for comparison with nu- 
clear collision data. The hard component from p-p spec- 
tra determines the initial parton spectrum parameters 
for nuclear collisions. Fig. [5] (left panel) shows a surface 



plot of the integrand of Eq. ^ — Dpp{y,ymax) 



dcrdi_ 



/cut(2/maa;) = {tanh[(y„aj; - j/c«t)/?cut] + l}/2 (6) incorporating FFs based on the LEP parametrization 
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plus the FF cutoff inferred from p-p collisions, p-p FF 
distributions are bounded below by Umin ~ 1-5 (pt ~ 0.3 
GeV/c). The plot z axis is logarithmic to show structure 
over the entire distribution support. 




FIG. 8: Left: (Color online) Argument of the pQCD folding 
integral on {y,ymax) based on in-vacuum p-p FFs. Right; 
Fragment distribution HNN-vac (integral on ymax) obtained 
from in-vacuum p-p FFs (solid curve) compared to the 
Gaussian-plus-tail model of the p-p hard component (dash- 
dotted curve) and the measured hard component from NSD 
p-p collisions (solid points) [3. The dotted curve is discussed 
in Sec. En] 

Fig. \E\ (right panel) shows the corresponding i?NN-vac 
FD (integration of the left panel over ymax) as the solid 
curve. The mode of the FD is ~ 1 GeV/c. The dash- 
dotted curve is a Gaussian-plus-tail model function, and 
the solid points are hard-component data from p-p colli- 
sions (l^ . The comparison determines parton spectrum 
parameters ycut = 3.75 {Ecut ~ 3 GeV) and exponent 
nqcD = 7.5. The data are well-described by the pQCD 
folding integral. 

The FD in the right panel represents the minimum-bias 
ensemble of jets that fall within the detector 77 accep- 
tance. The FDs are plotted as (1/y) cPrih/dy drj for com- 
parison to spectrum hard components plotted on trans- 
verse rapidity yt- The region above 2 GeV/c dominated 
by the pQCD power law is the conventional focus for 
study of parton fragmentation. Ironically, this analysis 
reveals that the physically most significant fragmentation 
structure and evolution hes below 2 GeV/c, the region 
conventionally assigned to hydro phenomena Q . 

Fig. [5] (left panel) shows the argument of the folding 
integral incorporating unmodified FFs from e-e collisions. 
The main difference is the extension down to ymin ^ 0.35 
{pt ~ 0.05 GeV/c). Fig. [5] (right panel) shows the corre- 
sponding FD (solid curve) . The parton spectrum param- 
eters determined by the p-p hard component are retained. 
The solid curve is the "correct answer" for an FD describ- 
ing inclusive hadrons from inclusive partons produced by 
free parton scattering from p-p collisions, which is not 
observed in real nuclear collisions (cf. Sec. IX Bp . The 
dash-dotted curve represents the hard-component model 
inferred from p-p collisions. The FD from e-e FFs lies 
well above the measured p-p hard component for hadron 
Pt < 2 GeV/c {yt < 3.3), and the mode is reduced to 




FIG. 9: Left: (Color online) Argument of the pQCD folding 
integral on {y,ymax) based on in-vacuum e'*'-e~ FFs. Right: 
Fragment distribution -ffeo-vac (integral on ymax) obtained 
from in-vacuum e^-e~ FFs (solid curve) compared to the 
Gaussian-plus-tail model of the p-p hard component (dash- 
dotted curve) [1^ . 



~ 0.5 GeV/c. The "correct" e-e FD strongly disagrees 
with the most relevant part of the p-p pt spectrum — 
the hard component. Despite strong disagreement the 
e-e FD is the correct reference for nuclear collisions, as 
shown below. 

The p-p spectrum hard component determines the par- 
ton spectrum cutoff energy, the only adjustable parame- 
ter in the folding integral since the shapes of the p-p FFs 
are defined by independent FF data. The p-p FF mul- 
tiphcity systematics in Fig. [HI suggest that the effective 
parton (gluon) spectrum cutoff is determined (for charge- 
particle measurements) by the requirement to produce at 
least one, and therefore two, charged hadrons. The cut- 
off inferred from p-p FFs {ycut 3.75) then provides an 
upper limit for e-e FFs, since the latter have substan- 
tially larger mean multiplicities. The caveat is especially 
relevant when modeling medium-modified FDs. 



D. NLO "fragmentation functions" 

NLO FDs (sometimes termed "fragmentation func- 
tions") [3^ are typically compared to the full p-p pt 
spectrum including the soft as well as hard compo- 
nents [37|- Theoretical representations of e-e FFs are 
currently less accurate for small fragment momenta (e.g.. 
Fig. 20 of 15]). As noted, inclusion of accurate in- 
vacuum e-e FFs in the pQCD folding integral leads to 
a large excess over the p-p hard component at smaller yt, 
which may have important physical significance (as dis- 
cussed below). The full soft-t-hard p-p spectrum greatly 
exceeds the FD at smaller yt- Comparison of NLO FDs 
with the full p-p spectrum can then suggest agreement 
with data which is misleading. Instead, they should be 
compared directly to the p-p spectrum hard component 
as in Fig. [5] (right panel), which then reveals physically 
important differences. 
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VII. THE PARTON SPECTRUM 

The pt spectrum for partons scattered from N-N col- 
lisions can be approximated by (\/pt) dadijet/dpt = 
Api "'J'^" above some cutoff Pt,cut- The three spectrum 
constants {A, uqcd, Pt,cut) can be inferred from nuclear 
collision data. Parton spectrum information comes from 
event-wise jet reconstruction and from the single-particle 
spectrum hard component. Jet reconstruction provides 
only a part of the differential spectrum. A and uqcd 
can be inferred from the larger-pt region. A pt-spectrum 
hard component combined with an FF ensemble can de- 
termine the parton spectrum cutoff explicitly from the 
structure near its mode. The parton (dijet) total cross 
section is then determined with improved accuracy. 



A. Parton spectrum from data 

When integrated over fragment rapidity y the folding 
integral in Eq. ^ can be expressed as the product of a 
weighted-mean dijet multiplicity and the integrated di- 
jet cross section. Since the parton spectrum is sharply 
peaked near 3 GeV the mean dijet multiplicity (Fig. [Al- 
right panel) is close to the value at that energy. Factor 
e is introduced to represent the ID dijet efficiency: the 
average fraction of a jet that falls inside 77 acceptance bin 
Sij (some jets overlap the boundary) and the fraction of 
partners in a back-to-back jet pair that also fall inside Srj. 
Both fractions depend on Srj, effective An interval A77 and 
the number N^ijet of jet pairs in a nuclear collision. 

The 7/t-integrated two-component model of mean 
hadron yields for NSD p-p collisions is 

drich _ dns_ ^ dnn 
drj drj drj 

The jet cross section can be inferred from measurements 
of spectrum hard component H = (1/y) d^nh/dy di]. In- 
tegrating Eq. ([3]) over y and expressing the remaining 
Vrnax integral as fidtjet (^dijet gives 

^ - |^^j.(^,,A,)n.,-. (9) 
d-q ycTNSD i^rj J 

for NSD p-p colhsions. The measured hard/soft ratio 
for NSD mean hch ~ 1-25 is nt/us ^ 0.008 Since 
dng/drj ^ 2.5 for NSD p-p collisions duh/drj ~ 0.02. 
Given gnsd ^ 36 mb, p-p fidijet ~ 3 and dijet fraction 
e ~ 0.45 (cf. Fig.[Tni) included in \r]\ < 0.5 {6r] = 1) then 
adijet/^V = 0.5 ± 0.12 mb. Assuming Stj^tt = ^ 5 
gives a total cross section adijet = 2.5 ± 0.6 mb. Uncer- 
tainty estimates are discussed in Sec. IXIII Equivalently, 
the probability of a minijet within \r]\ < 0.5 in NSD p-p 
collisions [expression within curly brackets in Eq. ([9])] is 
0.5 mb/36 mb ~ 0.014±0.003, consistent with [ij based 
on eudijet ~ 2.5 ± 1.0 and dnh/drj afichdrich/drj ^ 
0.03 ±0.01. 



The correspondence between a dijet cross section and 
average jet fragment yield in an 77 acceptance is non- 
trivial. The general result for A-A collisions depends 
on mean dijet number Ndijet = nunaryf^dijetlcrNSD- 
pQCD calculations produce differential cross section 
d^adijet/dpt dyi dy2, with dadijet/dpt or dadijet/dymax as 
a straightforward ID projection. To obtain a 2D projec- 
tion onto ymax and y or rj the single integral over y must 
accommodate the integer number of dijets in a collision. 

Fig. (TU] (left panel) shows calculated mean jet mul- 
tiplicities Njet for occupied bins of width Srj within 47r 
acceptance Arj. Nj^t varies as 2 Ndijet Srj/Arj toward the 
right. Toward the left where Ndijet ^ 1 the mean jet 
number in an occupied bin is 1/(1 — 0.5 Srj/ Arj). 




FIG. 10: Left: Mean jet number Njet in an occupied bin Sr] 
within acceptance A?; for dijet number Ndijet and assuming 
Poisson statistics. Right: Fraction e of dijet fragment yield in 
occupied bin Si] within interval A-q for dijet number Ndijet- 

In Fig. [To] (right panel) e{Ndijet,Sr], A77) represents the 
fractional yield of hadrons per dijet into acceptance Srj. 
For a single dijet in Arj ~ 5 (e.g., some p-p collisions) 
and Sri = 1, e(l,l,5) - 0.8/(2 - Sr]/Ar]) ~ 0.45. For 
Ndijet ~ 50 (e.g., semi-central Au-Au collisions) and 
Srj = 2, e(50, 2, 5) ~ 0.9, a two-fold increase in fractional 
hadron yield per dijet. 

Figure [TT] (left panel) shows the calculated p-p FD for 
spectrum cutoff j/cut = 3.75 (solid curve) and for j/cut = 
3.65, 3.85 corresponding to 10% shifts in the cut energy 
about 3 GeV (dotted curves) illustrating the precision of 
the cutoff determined from data: p-p data determine the 
cutoff to better than 5%. 



B. Parton spectrum from theory 

Figure [TT] (right panel) compares the spectrum defined 
by this analysis (solid curve, and note the factor 3) with 
theory and 200 GeV UAl data. As noted, the spec- 
trum from this analysis integrates to 2.5 ± 0.6 mb with 
well-defined cutoff ~ 3 GeV. The KLL parametrization 
600/p^ mb/(GeV/c) [H] (dashed fine) integrates to 2.2 
mb above 3 GeV/c. 

The UAl spectrum (points) is comparable at larger pt 
with the spectrum from this analysis (multiplied by 3) in 
an interval where the jet-finding efficiency should be good 
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FIG. 11: Left: HNN-vac for a nominal 3 GeV spectrum cutoff 
(solid curve) and ycut shifted by 0.1 {Ecut changed by 10%) 
in either direction (dotted curves) compared to p-p reference 
Hgg (dash-dotted curve) and NSD p-p data (solid points), 
illustrating the precision of the parton spectrum cutoff deter- 
mination. Right: The parton spectrum defined in this analy- 
sis with nqcD ~ 7.5 integrating to ~ 2.5 mb multiplied by 3 
(solid curve), KLL parametrized spectrum [18|] integrating to 
2.2 mb above 3 GeV/c (dashed curve) and UAI measured jet 
differential cross section integrating to ~ 4 mb (points) [19l |. 



and underlying-event contributions are relatively small. 
The slopes are the same but the amplitudes are different. 
At smaller pt the UAI data fall below the solid curve in a 
region where jet-finding efficiency rnight be reduced. The 
UAI spectrum integral is 4 mb [l9|. 

Several theoretical calculations of the jet total cross 
section stimulated by the UAI minijet data were intended 
to extend pQCD comparisons with data down to '--^ 3 
GeV. In [3^ a "two-component" spectrum model for p-p 
collisions was discussed, and it was estimated that be- 
cause of the underlying event the Et ^ 5 GeV observed 
by UAl corresponds to 3-4 GeV/c parton momentum. 
The total cross section obtained from a pQCD calcula- 
tion was 2-2.5 mb for a 3 GeV spectrum cutoff. 

A similar calculation in [2^ obtained 3 mb for a jet 
energy threshold of 3 GeV. However, the assertion that 
there is no actual threshold for parton — > hadron is con- 
tradicted for charged hadrons by the present study. The 
mean number of jets in NSD p-p collisions was given as 
'^cFdijet / cTineiastic^N SD ~ 0.1-0.2, which compares with 
2 X 2.5 ± 0.6 mb/ 36 mb = 0.14 ± 0.03 from this analysis. 

An extensive study of effects on the minijet cross sec- 
tion from energy scaling of structure functions was de- 
scribed in [l^]. The goal was to apply pQCD to low-pt jet 
physics — "the so-called minijet regime." Estimation of 
the underlying-event contribution to the UAl Et cluster 
finder lead to the conclusion that a 5 GeV Et calorimeter 
cluster (minijet) translates to a 3 GeV parton. Almost 
all partons in the spectrum (i.e., with energies ~ 3 GeV) 
are gluons. The minijet cross section was estimated to 
be 2-3 mb at 200 GeV, with spectrum cutoff near 3 GeV. 

Based on comparisons of spectrum hard components 
with calculated FDs and with pQCD theory there can 
be considerable confidence that the measured spectrum 
hard component from p-p collisions is a parton fragment 
distribution, and the corresponding minimum-bias par- 



ton spectrum is well defined. Given that baseline we 
now consider parton "energy loss" in A- A collisions. 



VIII. PARTON "ENERGY LOSS" 

In conventional pQCD descriptions of parton energy 
loss the leading parton is said to lose energy by gluon 
bremsstrahlun g du ring passage through a (possibly col- 
ored) medium [401 ■ In the context of this analysis some 
questions emerge: 1) Do complete FDs in more-central 
A- A collisions actually indicate energy loss (i.e., do FFs 
integrate to reduced parton energy)? 2) Do jet angular 
correlations manifest structure changes consistent with 
leading-parton random multiple scattering (e.g., broad- 
ening symmetric about the jet axis)? 3) If the answer to 
1) is "yes" how and where is the lost energy manifested 
in the medium? 4) If the answer to 1) or 2) is "no" what 
is the relevance or proof of an independent medium? To 
address those questions we can incorporate energy-loss 
models into calculated FDs from e-e and p-p collisions 
and compare with hard components (single-particle spec- 
tra and correlations) measured in Au-Au collisions. 



A. Negative boost of the p-p hard component 

A simple algebraic model of energy loss in FDs can 
be obtained by shifting p-p or N-N hard-component 
model Hqq down on rapidity (negative boost Aj/t ~ 



The model manifestly does not con- 
serve energy (energy /S.E is lost from the fragment sys- 
tem). Ratio rAA lS| is then modeled as hi{rAA) by 



In 



Hociyt) 



-Ayt , 10 

dyt 



where Hgg is the Gaussian plus power-law tail inferred 
from p-p collisions (denoted by n/i in [14*1 ) . The Hgg 
reference is included in plots of rAA below. The negative- 
boost model is imperfect because it does not respect the 
details of QCD splitting, and the lower limit of the FD 
experiences the same rapidity shift as the leading parton, 
which is probably not physical (cf . the next subsection) . 
Its recommendation is algebraic simplicity. 



B. "Medium-modified" FFs 

A better-justified model (BW) of medium-induced 
modification to QCD fragmentation is described in [29j . 
Parton "energy loss" is modeled, but the parton energy 
is conserved within the modified FF. In contrast to spe- 
cial treatment of the leading parton (bremsstrahlung) 
with "broadening and softening" of the FF, all sublead- 
ing splittings are treated equally by BW — momentum is 
conserved at all stages of the cascade. 

The BW model is applied to MLLA descriptions of two 
experimental FFs (TASSO 14 GeV [Ml and OPAL 200 



GeV [ill). As noted in [T5| MLLA FF parametrizations 
fail for small and large fragment momenta. The bottom 
20% of the FF, where the most interesting conclusions 
of this analysis emerge, is typically missing. Discrepan- 
cies are typically of the same magnitude as the medium 
effects observed in this analysis using an accurate FF rep- 
resentation. The statement "the MLLA can serve as a 
baseline in searching for medium effects" is not justified. 

Figure [TJlfleft panel) illustrates the BW model (cf. 
Fig. 1 of |2i]). In- vacuum FFs for Q = U and 200 
GeV derived from the beta parametrization are shown 
as the dashed and solid curves respectively [l^. Data 
from TASSO 14 GeV (solid points) are shown for com- 
parison, duplicating part of Fig. [3] (left panel). The 
beta distribution represents all FF data to the statistical 
limits. The practical consequence of the BW "energy- 
loss" mechanism is a momentum-conserving rescaling of 
FFs on momentum fraction Xp or logarithmic variable 

= hi{l/xp). Density reductions at larger fragment mo- 
menta (smaller ^p) are balanced by much larger increases 
at smaller momenta. The large changes correspond to an 
inferred leading-parton fractional "energy loss" of 25%. 




= ln(Pje/P) y = ln[(E + p)/mj 

FIG. 12: Left: parametrized fragmentation functions from 
e^-e~ collisions at two energies plotted on for in-vacuum 
FFs (dashed and solid curves) [1^ s^nd for in-medium modifi- 
cation by rescaled splitting functions (dash-dotted and dotted 
curves) p9l | compared to data for 14 GeV (solid points) [sH - 
Right: Curves in the left panel replotted on rapidity y. 

The BW model modifies the splitting process by rescal- 
ing the momentum fraction consistently at all stages of 
the cascade. The beta parametrization of FFs [l^ has 
only two parameters {p, q) . p represents the effect of 
quantum coherence or gluon saturation which effectively 
terminates the cascade at hadron formation (LPHD). q 
represents the pQCD splitting cascade itself. Empiri- 
cally, we observe that increasing parameter g by 0(1) in- 
crement Aq accurately duplicates the BW rescaling pro- 
cess. 

The dotted and dash-dotted curves in Fig. [T^] corre- 
spond to q q + 1.15. The curves (left panel) match 
Fig. 1 of [29| in the interval for which the MLLA approx- 
imates FF data. We can thus implement the BW energy- 
loss model simply by changing q in the beta parametriza- 
tion to achieve an accurate and complete representation 
of measured FFs. Figure d] (right panel) uses the beta 
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parametrization to demonstrate the multiplicity increase 
(open circles) corresponding to the BW "energy-loss" 
(medium modification) prescription. 

Figure[l2] (right panel) shows the same system on frag- 
ment rapidity y. The relation is given by — lii{pjet/p) 
= h\{2 pjet/m-^) - h\{2p/m^) ~ y^ax - y, with energy 
scale Q = 2pjet- Beta FFs on don't extend to infinity 
because FFs on y are bounded below by ymin- In e-e 
collisions ymin ~ 0.35 corresponds to p ^ The 
maximum ^p value is thus ip.max ^ ln(Q/m^) = y„iax 
(~ T in [13). 

The BW "energy-loss" model coupled with the FF beta 
parametrization from fisj provides an accurate algebraic 
model of FF "medium modification" valid for all rele- 
vant energies and momenta and directly related to pQCD 
principles. The system can generate medium-modified 
FDs for comparison with A-A hard-component evolution. 



C. "Energy loss" and FDs 

Figure [13] (left panel) shows the e-e FF ensemble with 
BW modification as described in the previous subsec- 
tion. Energy-loss parameter is the change in beta- 
distribution model parameter q which emulates the BW 
energy-loss method. The value Ag — 1.15 (for 0-12% cen- 
tral Au-Au collisions) is determined by hard-component 
data above pt ^ 4: GeV/c, where "suppression" in ratio 
rAA is approximately constant. For this initial survey Aq 
is assumed to be independent of parton energy (ymax)- 




FIG. 13: Left: (Color online) Argument of the pQCD folding 
integral on {y,ymax) based on in-medium e+-e~ FFs. Right: 
Fragment distribution Hcc-mcd (integral on ymax) obtained 
from in-medium e^-e~ FFs (solid curve) compared to the in- 
vacuum FD (dotted curve) and the Gaussian-plus-tail model 
of the p-p hard component (dash-dotted curve) [T3 ]. 

Figure [T3] (right panel) shows i/co-mod (solid curve) , 
the FD obtained by inserting e-e in-medium FFs from 
the left panel into Eq. ([3]) and integrating over parton 
rapidity ymax- The dotted curve is the i/cc-vac reference 
from in-vacuum e-e FFs. The dash-dotted curve is again 
the Gaussian-plus-tail p-p hard component Hqq for ref- 
erence. The mode of iJoo-mod is ~ 0.3 GeV/c. 

Figure [HI shows results for p-p FFs. The major dif- 
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^max 



The soft component is interpreted as longitudinal pro- 
jectile fragmentation [via inelastic N-N scattering) ap- 
proximately independent of A-A centrality. The hard 
component is interpreted (as in p-p collisions) as the 
FD from minimum-bias parton scattering and fragmen- 
tation into some hadron angular acceptance. Because a 
single soft-component model function is subtracted from 
all centralities any systematic error in the subtraction is 
common to all hard components. Relative variations with 
centrality are then unique to hard-component structure. 



FIG. 14: Left: (Color online) Argument of the pQCD folding 
integral on {y,ymax) based on in-medium p-p FFs. Right: 
Fragment distribution -ffNN-med (integral on i/max) obtained 
from in-medium p-p FFs (solid curve) compared to the in- 
vacuum FD (dotted curve) and the Gaussian-plus-tail model 
of the p-p hard component (dash-dotted curve) [ij . 



ference between p-p and e-e FDs appears below pt ^ 2 
GeV/c (yt ~ 3.3). Conventional comparisons with theory 
(e.g., data vs NLO FDs) typically do not extend below 2 
GeV/c [13. 



B. Hard-component spectra and ratios 

Figure [T51f left panel) shows hard-component evolution 
with centrality for pions from 200 GeV Au-Au collisions 
(five centrality classes). The spectrum data are in the 
form of a 3D density on (7/4,77,0). The thin dotted ref- 
erence curves are obtained from Eq. pip by replacing 
Haa with reference Hnn (model function Hgg)- The 
points are the hard component from 200 GeV NSD p- 
p collisions ll^]. The main features are the suppression 
at larger yt intensively studied at RHIC and the much 
larger enhancement at smaller yt first described in 



IX. A-A TWO-COMPONENT SPECTRA 

This analysis describes spectrum hard components by 
folding a parton power-law spectrum with parametrized 
FF ensembles. The hard component from p-p collisions 
constrains the parton spectrum. Hard components have 
also been extracted from Au-Au spectra for five centrali- 
ties [1| . A sharp transition in hard-component properties 
is observed as for minijet angular correlations [l3|. We 
use measured Au-Au hard components to study the evo- 
lution of parton fragmentation (and possibly the underly- 
ing parton spectrum) with centrality in nuclear collisions. 



A. A-A two-component model 

The algebraic form of the two-component model of per- 
participant-pair A-A spectra is 



2 1 drich 
ripart yt dyt 



= SNNiVt) V HAAiuuv) 



(11) 



SNNiUt) + vrAA{vu ^) HNNiVt), 



where Snn (~ Spp) is the soft component and Haa is the 
hard component (with reference Hj^n ^ Hpp) integrat- 
ing respectively to multiplicities ris and rih in one unit 
of pseudorapidity r] [1, Ratio taa — Haa/ Hjyjy is a 
refinement of nuclear modification factor Raa- Central- 
ity measure f = 2nunary / riparUcipant estimates a mean 
participant-nucleon path length in the Glauber model. 
Model functions Sxx and Hxx are normalized to be com- 
patible with the total spectrum density in a given context 
(i.e., whether the density is ID, 2D or 3D). 
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FIG. 15: (Color online) Left: Spectrum hard components 
from five centralities of 200 GeV Au-Au collisions (bold curves 
of several types) 01 compared to a two-component reference 
system (thin dotted curves). Hard-component data from p- 
p collisions (solid points) [ij] are included for comparison. 
Right: Hard-component ratio taa for the data in the left 
panel (bold curves of several types and dots) compared to a 
simple energy-loss model (black dash-dotted curve). 

Figure [T5l (right panel) shows corresponding ratio rAA 
based on hard-component reference Hnn set equal to 
Gaussian model Hqg = UhHo from [14]. Evolution of 
suppressions and enhancements is more clearly visible. 
The p-p data and the most peripheral Au-Au data agree 
with the N-N reference {vaa — 1) above yt — 2.5 but 
deviate significantly from Hgg below that point. 

As noted in [8] there is a sharp transition in the cen- 
trality trend for both suppression and enhancement, also 
seen in 200 GeV Au-Au angular correlation studies p^ . 
The direct correspondence between the centrality trend 
at pt ~ 10 GeV/c and that at pt ~ 0.3 GeV/c seems cu- 
rious in the conventional RHIC context (hydro -I- high-pt 
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jets). However, if the entire hard component is inter- 
preted as an FD the correlation is seen to be inevitable. 

A major consequence of this analysis is the realization 
that HjyN from p-p collisions is not the correct FD ref- 
erence for nuclear collisions. The p-p hard component is 
itself strongly modified relative to the correct reference. 
An alternative to vaa is required for differential study. 



X. CENTRALITY EVOLUTION OF THE FD 

The insensitivity of nuclear modification factor Raa 
to most energy-loss details and the superiority of hard- 
component (FD) ratio rAA = Haa/H^n were demon- 
strated in [g]. In this section calculated FDs are com- 
pared directly with measured A-A spectrum hard com- 
ponents using generalized hard-component ratio r^x to 
determine centrality evolution. The notation adopted 
is FD -ffxx, where xx = pp (data), A A (data), GG 
(Gaussian-plus-tail model) and NN-vac, NN-med, ee-vac, 
ee-med. The last four, introduced in the previous section, 
are obtained from folding integrals. We then define ratios 
TAA Haa/Hgq (used in [1] and Fig.[Ti]- right panel). 

Toe = i?cc-mcd/-ffoc-vac and TNN = ^NN-mcd/^NN-vac . 



A. FD ratios compared to central Au-Au collisions 



Figure [TBI (right panel) introduces a novel concept. In- 
stead of comparing the calculated in-medium FD for N-N 
(averaged within A-A collisions) with the in-vacuum FD 
for N-N or similarly comparing e-e with e-e as in the left 
panel, the in-medium FD for e-e is compared with the 
in-vacuum FD for N-N by defining ratio 



NN-vac 



(12) 



Calculated r^N describes the measured taa well over the 
entire fragment momentum range. The question then 
arises how to interpret the result. 



B. Revised FD reference for nuclear collisions 

Fig. [16] (right panel) implies that the FD reference for 
nuclear collisions should be reconsidered. The hard com- 
ponent for central Au-Au collisions appears to be well- 
described by i?cc-mcd with the BW "energy loss" mecha- 
nism applied to e-e FFs. The p-p hard component devi- 
ates strongly from H^c-vac in Fig- El (right panel) and may 
be strongly suppressed at smaller yt- The combination 
suggests that p-p hard component Hnn is not the cor- 
rect reference for A-A collisions as assumed implicitly in 
defining conventional ratio Raa- The proper in-vacuum 
reference for nuclear collisions is actually i?cc-vac- 



Figure [16] (left panel) shows calculated FD ratios rpf^ 
(dashed curve, p-p FFs) and r^e (dash-dotted curve, e-e 
FFs). The solid curve is the measured vaa from central 
(0-12%) Au-Au collisions at 200 GeV Ij. Aq - 1.15 
for Hee-med and H Ai N -med was adjusted to obtain the 
correct large-j/t suppression for that centrality. The ref- 
erence for rAA is hard-component model function Hqg- 
The dotted curve is the ratio reference obtained by shift- 
ing Hqg on yt by Aj/t '-^ —0.26 (negative boost). As 
noted in Q the simple negative-boost model does not 
describe the Au-Au data. But the e-e and N-N ratios 
also do not describe the data. 
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FIG. 16: Left: Calculated hard-component ratios r^x com- 
pared to measured rAA for 0-12% central 200 GeV Au-Au 
collisions (bold solid curve). Right: Same as the left panel 
except for newly-defined ratio reN ~ J^cc-mcd/^^^NN-vac which 
compares well with the Au-Au data. 
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FIG. 17: (Color online) Left: Hard-component ratios r^x 
based on the calculated reference Hee-vac determined with in- 
vacuum e^-e~ FFs. The data are for Hpp from p-p collisions 
(solid points) and Haa from Au-Au collisions with v < 2.5 
(solid curve), ree includes the calculated //ee-med from in- 
medium e"'"-e~ FFs (dash-dotted curve). rNN includes the 
calculated i/NN-vac from in-vacuum p-p FFs (dashed curve). 
Right: Same as left panel but for Haa from Au-Au collisions 
with u > 2.5. The dotted curves are discussed in the text. 

Figure [T7| (left panel) shows FD ratios redefined in 
terms of the ee-vac reference: Hpp (p-p data - points), 
-ffAA (peripheral Au-Au data ~ solid curve) and calcu- 
lated i/oo-mcd (dash-dotted curve) and i/NN-vac (dashed 
curve) all divided by reference i/co-vac- The strong sup- 
pression of p-p and peripheral Au-Au data apparent at 
smaller yt results from the cutoff of p-p FFs noted above. 

Figure [T71 (right panel) shows measured Haa/ Hee-vnc 
for more-central Au-Au collisions (solid curves) above a 
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transition point on centrality at v ~ 2.5. Centrality mea- 
sure V = 2ni,in/npart IS the mean participant-nucleon 
path length in number of N-N coUisions. For the Au-Au 
cohisions in Fig. [T7| values for the five centralities are 
1.93, 2.83, 3.92, 4.87, 5.5, where ~ 1.25 is N-N colli- 
sions and 1/ ^ 6 is 5 = Au-Au collisions [1, li^l . From 
= 1.98 to ly — 2.83 there is a dramatic change in the hard 
component. At the transition point j/ ~ 2.5 Upart = 40 
(out of 382) and nun = 50 (out of 1136) 



TABLE I: Parameters for _ffxx vs centrality 
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centrality 


Aq 


yo = ?y 


ycut 


Ecnt (GeV) 


XX 


reference 


0.0 


0.0 


3.75 


3.0 


ee-vac 


1 


0.0 


1.5 


3.75 


3.0 


NN-vac 


2 


0.7 


1.25 


3.7 


2.85 




5 


1.15 


1.25 


3.65 


2.7 




reference 


1.15 


0.0 


3.65 


2.7 


ee-med 



Table [J shows parameters for calculated FDs H^^ vs 
Au-Au centrality plotted as ratios to common reference 
iJee-vac in Fig. [T7] (right panel). Aq is determined in all 
cases by the region above yt — 4: (conventional "sup- 
pression"). Ucut or Ecut is determined by the slope in 
the intermediate region near yt = 3. FF cutoff param- 
eters yo — are determined by structure to the left 
of yt = 2.5. The first line describes reference i?cc-vac 
{rxx = 1) with cutoff energy Ecut — 3 GeV. The second 
line describes H^N-vac (dashed curve). The last line de- 
scribes the limiting case of -ffco-mcd for central Au-Au col- 
lisions (dash-dotted curve) with Ecut — 2.7 determined 
by the data near yt = 3. The lighter dash-dotted curve 
is the same with Ecut — 3 GeV for comparison. Dotted 
curves 2 and 5 are hybrid versions -ffxx-mod with yo — 
adjusted to accommodate the data to the left of yt = 2.5. 

The results can be interpreted as follows. With in- 
creasing centrality the splitting cascade is modified {Aq 
increases from zero), suppression of FFs at smaller y is re- 
duced (j/o = rnove toward zero) and parton spectrum 
cutoff ycut {Ecut) is also reduced, increasing adijet and 
the total minijet yield by 50% as in Fig. [7] (left panel). 
Instead of invoking r^x ratios the measured spectrum 
hard components can be compared directly with calcu- 
lated FDs to reveal fragmentation evolution. 



C. Fragmentation evolution 

Figure [18] shows spectrum hard components Haa 
(solid curves) for five centralities from 200 GeV Au-Au 
collisions |8|. This format shows Haa and related FD 
curves, whereas Fig. [15] (left panel) shows v Haa includ- 
ing participant path length v. Since vupart/^ = ntmary 
the hard components of un-normalized yt spectra (2D 
densities) scale proportional to ribinary as expected for 
parton scattering and fragmentation in A- A collisions. 
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FIG. 18: (Color online) Measured spectrum hard components 
Haa for five centralities from 200 GeV Au-Au collisions (bold 
curves of several colors) and 200 GeV NSD p-p collisions 
(solid points) compared to calculated FDs for several condi- 
tions (vacuum, medium, e^-e~, p-p). The hatched region at 
upper left esimates the uncertainty due to the 5'jvjv subtrac- 
tion common to all centralities. The hatched region at lower 
right denotes the interval conventionally allotted to pQCD. 



The points are hard-component data from 200 GeV 
NSD p-p collisions 14]. The dash-dotted curve is the 
standard Gaussian -|-tail model function Hqq. FDs from 
the previous section are also shown. The dashed curve is 
^^NN-vac, and the upper (bold) dotted curve is ^^oc-mod 
with Aq — 1.15, which nominally corresponds to the 
most-central Au-Au curve (0-12%). The parton spec- 
trum cutoff for Hee-med has becu reduced from 3 GeV 
{Vrnax = 3.75) to 2.7 GeV {ymax = 3.65) to match the 
central Au-Au hard component near yt = 3. The two 
thinner dotted curves labeled 2 and 5 (Au-Au centrali- 
ties) are -ffoc-mod with cutoff parameters yo = S.y reduced 
to accommodate the data below yt = 2.5. 

The p-p and peripheral Au-Au data are consistent with 
^^NN-vac by construction. Above a transition point on 
centrality {v ^ 2.5) Haa transitions from ifNN-vac to- 
ward i?oc-mcd over most of the yt range, with residual 
deviations confined to smaller yt- For more-central colli- 
sions agreement of Haa with the limiting -ffoc-mcd curve 
extends toward the limits of accepted yt- 



D. Restoration of the FF base in A-A collisions 

Below centrality i/ = 2.5 we observe strong fragment 
suppression at smaller yt relative to the -ffco-vac reference. 
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Above that point there is asymptotic approach to i?cG-mcd 
(bold dotted curve in Fig. [T5)) . with strong enhancement 
at smaller yt and the expected suppression at larger yt 
characterized as "jet quenching". The FD modification 
in central Au-Au collisions is approximately consistent 
with the BW in- medium modification of e-e FFs [2^ . 
Those trends suggest the following scenario: 

• The bases of fragmentation functions (jets) in p-p 
and peripheral A-A collisions are missing from re- 
constructed jets and the spectrum hard component 
compared to in- vacuum e-e FFs 

• Above a transition point on A-A centrality frag- 
mentation changes dramatically 

• The bases of A-A FFs are partially restored to com- 
patibility with e-e FFs 

• A-A FFs are modified in a manner compatible with 
a pQCD description of FF medium modification 

• Most of the underlying parton spectrum does not 
change with Au-Au centrality; no scattered partons 
are lost to the final state (thermalized) 

• However, the low-energy cutoff — 3 GeV for N-N 
collisions — falls to 2.7 GeV for central Au-Au col- 
lisions, increasing the minijet cross section by 50% 

The transition of the Au-Au spectrum hard component 
near u ~ 2.5 corresponds to a similar sharp transition 
observed in minijet angular correlations at 200 GeV p^ . 
again consistent with equivalence of the spectrum hard 
component and minimum-bias parton fragmentation. 

XI. RELATED ASPECTS OF A-A COLLISIONS 

Given a better understanding of scattered-parton spec- 
tra and fragmentation we consider the consequences in 
A-A collisions for parton distributions in the initial state 
and hadron distributions in the final state. 



A. pQCD parton spectra and A-A initial conditions 

Parton spectra determined ah initio from pQCD are 
used to estimate the initial conditions for heavy ion colli- 
sions at RHIC. In [i^ a jet spectrum was obtained from 
a parton spectrum defined in terms of a pQCD parton 
differential cross section and nucleon structure functions 

dn^ = KT{0)^^, (13) 
apt dpt 

where K = 2 ("higher-order contributions") and T(0) 
is the "nuclear geometrical factor" = dA^ /8tt R\, with 
Ra = I.IA^/^ and T(0) = 34 mb'^ for central Au-Au 
collisions. The pQCD spectrum was determined as in this 



analysis for central Au-Au collisions at 200 GeV where 
nbinary — 1136 and unsd — 36.5 mb. T(0) can then 
be compared to nunaryhyNSD — 31 mb^^. Extracting 
dujet/dpt data from ^431 (Fig. 1, 2-10 GeV) and rear- 
ranging Eq. to 

dcJdijet _ Pt drijet 
dy max 

KT{0) dpt ' ' 

[with yraax = ln(2p(/m7r)] we recover the pQCD par- 
ton spectrum, plotted as the bold dotted curve in Fig. [7| 
(left panel) , corresponding closely (near the peaks) to the 
spectra inferred from FD data. 

Whereas FD data and this analysis imply that the 
parton spectrum terminates near 3 GeV with a 2.5-4 
mb total cross section, the spectrum in [43] was inte- 
grated down to 1 GeV to estimate a parton (minijet) 
density drijet /dy^ = 750 for central Au-Au collisions. 
The lower cutoff was justified by saturation-scale argu- 
ments [iil (cf. App. IXIII"G1) . Given in rapidity interval 
= 7, Ubinary 1136 and two jets per parton collision the 
implied dijet total cross section is about 85 mb (and cf. 
Fig. [7]- left panel), larger than the ioiaZ N-N cross section 
and more than 20 times the value we infer from hadron 
spectrum data. Based on a 4 mb jet cross section for cen- 
tral Au-Au collisions we expect 35±9 minijets in one unit 
of ?7, consistent with minijet correlation analysis ^10ll2l|. 

Arguing by analogy, a parton spectrum observed via 
charged hadrons should be terminated by the available 
density of hadronic final states, the same mechanism that 
terminates a splitting cascade in jet formation. An iso- 
lated parton scatter (e.g. in 200 GeV p-p collisions where 
at most one parton scatter occurs) should not proceed 
unless there is at least one hadronic final state available. 

It could be argued that in more-central Au-Au colli- 
sions the parton spectrum is substantially altered (cutoff 
extended to much lower energies) by the environment, 
that partons scatter and rescatter until hadrons finally 
emerge from a collective medium (cf. Sec. IXIII G[) . But 
that is not what we observe in direct comparisons of cal- 
culated FDs with measured spectrum hard components 
and correlations. We do see a modest decrease (10%) in 
the cutoff energy with corresponding 50% increase in the 
minijet cross section, as in Fig. [7] (left panel, dash-dotted 
curve). However, strong constraints from the hadron 
density of states apparently still apply to individual par- 
ton scatters. Rescattering of partons and hadrons or 
"constituent quark recombination" from a medium are 
contradicted by hard-component data. 



B. Charged-hadron and total-pt production 

Models for initial-state parton scattering and fragmen- 
tation should confront measured features of the final 
state. Beyond spectrum hard components this analy- 
sis can be compared with soft and hard components of 
integrated charged-hadron yields and total pt- The two- 
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component model for A-A particle production (with im- 
pact parameter b and participant path-length v) is 



2 drich dUg duh 



Upart drj 



drj drj 

V <Jdijet{b) 



(15) 



= 2.5-, 

(^NSD Ary 
= 2.5 + 0.02 p-p collisions 

« 2.5 + 1.1 6 = Au-Au collisions, 

defining dnh/drj (integral of Haa on yt) as the average 
hard component for a single N-N collision within an A- 
A collision. The soft-component multiplicity retains the 
p-p value 2.5. The jet cross section (and therefore minijet 
number) increases by about 50% with A-A centrality due 
to reduction of the parton spectrum cutoff, and the mean 
dijet multiplicity increases about 3-fold. In those A-A 
collisions with multiple jet pairs eaa ~ 2epp, effectively 
doubling the number of observed jets per N-N collision. 
A similar expression for the total-pt density is 



^part 



Pt 



^(.). + .^(.).(.) 



(16) 



= 0.88 GeV/c -I- 0.02 GeV/c p-p 

w 0.88 GeV/c + 0.5 GeV/c 6 = Au-Au, 

where Pt is the total pt in one unit of 77, {pt)s 0.35 
GeV/c is fixed, and {pt)h{b) is 1 GeV/c for p-p col- 
lisions but decreases with increasing A-A centrality to 
about 0.5 GeV/c due to medium modification of the FFs. 
Note that 1.38 GeV/c / 3.6 = 0.38 GeV/c, consistent 
with the {pt) centrality variation for Au-Au collisions [1]: 
{pt) increases with A-A centrality through a maximum, 
and then decreases due to the FF medium modification. 

In contrast, the hard-component contribution from (43l | 
predicting drijet/dy = 750 (parton spectrum cutoff at 1 
GeV) for central Au-Au collisions at 200 GeV should be 



^ - X — X 1 GeV/c 

"^part 3 191 



(17) 



2.5 GeV/c 5 = Au-Au coUisions, 

exceeding by five times the hard-component Pt observed 
in data. Measured charged-hadron minijets correspond- 
ing to a 4 mb jet cross section account for all hard- 
component Pt production. That reasoning closes the loop 
among correlations, single-particle spectra and pQGD. 



XII. SYSTEMATIC UNCERTAINTIES 

This analysis compares calculated FDs with measured 
spectrum hard components over a large kinematic range. 
The critical elements of the comparison are the parton 
spectrum, parametrized FFs from e-e and p-p collisions, 
an "energy loss" model applied to FFs, and spectrum 
hard components from p-p and Au-Au collisions. 



A. Parton spectrum 

In this analysis a jet cross section is obtained by com- 
paring a pQCD calculation with a minimum-bias frag- 
ment distribution rather than with a reconstructed-jet 
spectrum. The jet differential cross section on 77 and its 
relative error can be extracted from Eq. [5] 



^ dijet 

Ar] 



da. 



dijet 



<7NSD 



drih 



dr] 



e{dri,Ari)ndijet dr] 



(18) 



The parameter values with estimated systematic errors 
are (Jnsd = 36.5 ± 2.5 mb [H, dnn/dr] = 0.02 ± 0.0025 
([11] and this analysis), Udijet = 3.0 ±0.5 (GDF FFs [l3 
and systematic studies for this analysis) and t{5r] = 1) = 
(0.8±0.1){l-|-0.1±0.025}/2 (estimate based on the an- 
gular widths of observed minijets and the small curvature 
of the away-side ridge within |7y| < 1 [l3|)- The combined 
error for the differential cross section is then about 25%. 

Given the basic power-law form with cutoff, the par- 
ton spectrum is defined by three parameters — ycut (Ecut), 
nqcD and amplitude Ay^^^. The three are simply re- 
lated to adijet in Sec. IVI 51 uqcd and amplitude Ay^^^^ 
are determined by comparison to the p-p hard compo- 
nent up to 7 GeV/c, as in Fig. [HI fright panel), and the 
peripheral Au-Au hard component up to 10 GeV/c, as in 
Fig. [m Values uqcd = 7.5 ± 0.5, = (8 ± 1) x 10^ 

and ycut = 3.75 ± 0.05 {Ecut = 3.0 ± 0.15 GeV) are con- 
sistent with Udijet = 2.5 ± 0.6 mb. 

While the spectrum cutoff for p-p collisions is well- 
defined by data the effective cutoff for e-e FFs is less 
well-defined. The cutoff used for the in- vacuum e-e refer- 
ence FD is by definition the same as inferred from the 
p-p hard component. The cutoff used for centralities 
above the transition in Au-Au collisions was reduced to 
3.65 {Ecut ^ 2.7 GeV) to match the corresponding FDs 
to Au-Au hard components in the region near yt ^ 3 
{pt — 1.5 GeV/c). The uncertainty in the modified cutoff 
energy determined by central Au-Au data is 0.15 GeV [cf. 
Fig. [T7| (right panel) thick and thin dash-dotted curves] . 



B. Fragmentation functions 

The beta parametrizations in [l^ describe e+-e^ FF 
data to their error limits from 10 to 200 GeV energy 
scale (dijet energy) and for all fragment momenta from 
0.1 GeV/c to the parton momentum. Measured abso- 
lute fragment yields (dijet multiplicities) are also well 
described. The parametrization allows reliable extrapo- 
lation down to Q = 4-6 GeV {Ejet = 2-3 GeV). 

The p-p FF parametrization used in this analysis is 
less-well developed. Measured FFs plotted on rapidity 
for the first time in [l^ revealed substantial systematic 
deviations at small fragment momenta (suppression) rel- 
ative to e~'"-e~ FFs, although both are reported to be 
well-described by the MLLA [l^ . This analysis sug- 
gests that the deviations are due to physical differences, 
not detection or reconstruction inefficiencies. 
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Proper description of the suppression of p-p FFs is 
essential for the comparison between calculated FDs and 
spectrum hard components. E.g., the suppression a,t Q — 
6 GeV amounts to more than a factor 2 decrease in dijet 
multiplicity. Yet the p-p FFs so described lead to good 
agreement with theoretical jet cross sections. 

A two-parameter tanh cutoff function with parame- 
ters {yo,£,y) applied to the e^-e^ FF parametrization de- 
scribes CDF FFs well, especially for smaller dijet ener- 
gies. The FF cutoff influences the shape of calculated 
FDs near the spectrum mode, as does the parton spec- 
trum cutoff i/cut- The parameters could interact, for ex- 
ample in describing the p-p hard component. In Fig. [6] 
(right panel) the dotted curve shows the result of reduc- 
ing parameter z/q from 1.6 to 1.4 (which would strongly 
disagree with measured p-j5 FFs). The difference can 
be compared with the result of changing ycut by 0.1 in 
Fig. [m (left panel) . The shape changes are easily distin- 
guished at the few-percent level. 

C. Energy-loss model 

A simple model of FF medium modification was 
adopted for this study. It has the advantage that al- 
teration of one parameter (g) in the beta-distribution de- 
scription of e~^-e~ FFs reproduces the BW pQCD model. 
An accurate description of "medium modification" over 
the complete range of fragment momenta is thus possible. 

"Energy-loss" parameter Aq determined by spectrum 
suppression at larger pt produces large effects at smaller 
Pt which were unanticipated but describe Au-Au data 
well. Aq modifies beta distribution control parameters 
(p, q) derived for e+-e" FFs 'ISj. A value Aq = 1.15±0.2 
for central Au-Au collisions describes suppression of the 
hard component up to 10 GcV/c to the error limits of 
data. The resulting e-e FFs for 14 and 200 GeV agree 
with the BW result to a few percent over the range of 
validity of the MLLA FFs used in the BW study 

D. Spectrum hard components 

Spectrum hard components are obtained as differences 
between measured spectra and soft component Sxx in- 
ferred from the multiplicity (in p-p collisions p^ ) or cen- 
trality (in Au-Au collisions [1]) variation of a spectrum 
ensemble. Systematic uncertainties derive from the mea- 
sure spectra and from the subtracted soft component. 

The limit procedure used to obtain the soft component 
(a form of Taylor expansion) suppresses systematic un- 
certainties derived from data and the inferred structure 
of Sxx in the difference defined as the hard component. 
What error remains is common to all extracted hard com- 
ponents. Differential centrality variation of hard compo- 
nents is then relatively free of systematic error because of 
the double suppression of common-mode error. The ab- 
solute uncertainty due to the Sxx subtraction, common 



to all centralities, is estimated by the upper-left hatched 
region in Fig. [THl 



XIII. DISCUSSION 

This analysis provides direct comparison of the com- 
plete hadron-spectrum hard component from nuclear col- 
lisions with parton fragment distributions derived per 
pQCD from fragmentation functions and a common par- 
ton spectrum. Evolution of the spectrum hard compo- 
nent with Au-Au centrality reveals that the FF ensemble 
undergoes a transition from in-vacuum p-p FFs to in- 
medium e+-e^ FFs at a specific centrality. Details and 
implications of the analysis are now discussed. 



A. Hard component as fragment distribution 

Spectrum hard components previously extracted from 
p-p and Au-Au collisions are identified as minimum-bias 
fragment distributions or FDs. Calculated FDs are in 
turn generated with improved accuracy by folding beta 
distribution FF parametrizations with a power-law par- 
ton spectrum model. Comparison with the p-p hard com- 
ponent determines a spectrum cutoff and minimum-bias 
jet cross section consistent with pQCD and theoretical 
analyses of UAl Et clusters. 

Such comparisons confirm in detail that a pQCD de- 
scription of parton scattering and fragmentation is valid 
down to at least 3 GeV parton energy, the effective lower 
limit for partons fragmenting to charged hadrons. Fur- 
ther confirmation comes from correlation hard compo- 
nents which reveal unambiguous jet angular correlations 
down to small hadron (~ 0.1-0.3 GeV/c). 

These combinatoric methods have several advantages 
compared to the UAl event-wise Et analysis. The UAl 
analysis imposed a jet definition (cone jet finder) with 
attendant bias. A substantial background came from the 
underlying event (1-2 GeV within the jet cone). In con- 
trast, the hard component of particle spectra reveals the 
lower edge of the FD with relatively small bias (the sub- 
tracted soft component is defined by a limiting procedure 
unrelated to hard-component interpretation). The par- 
ton spectrum cutoff is thus accurately determined. 



B. New fragmentation phenomenology 

We observe that p-p FFs are systematically different 
from e-e FFs (e.g., OPAL/TASSO vs CDF), but the dif- 
ference is apparent only when FF data are plotted on 
rapidity y down to small fragment momenta as in p^ . 
The hard component of 200 GeV p-p collisions is well 
described by folding parametrized p-p FFs with a power- 
law parton spectrum terminating near 3 GeV. Since the 
p-p minimum-bias hard component is well described by 
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p-p FFs jet-reconstruction and particle-detection ineffi- 
ciencies seem not to cause the differences. 

We tlierefore conclude that p-p FFs are indeed sup- 
pressed at smaller fragment momenta relative to e-e FFs, 
implying that novel QCD physics emerges in p-p colli- 
sions, and the p-p hard component is not the proper FD 
reference for nuclear collisions. Use of the p-p FD (or 
worse, the entire p-p pt spectrum) as a reference appears 
to be misleading. The correct FD reference for all nu- 
clear collisions is obtained by folding in-vacuum e-e FFs 
with a power-law parton spectrum. A unique aspect of 
the present analysis (compared to conventional NLO FD 
calculations which follow the same approach) is the ac- 
curacy and extent of the e-e FF parametrization in [l^ . 

The centrality evolution of FDs is studied in this anal- 
ysis directly with spectrum hard components and in- 
directly with ratio measures. A new phenomenon has 
emerged: At a particular centrality there is a sharp tran- 
sition from the suppressed p-p FD to a medium-modified 
e-e FD with at least three-fold increase in dijet multi- 
plicity. The p-p FD observed in peripheral Au-Au be- 
comes in central Au-Au not an "energy-loss" p-p FD but 
a "medium- modified" e-e FD, and the effective cutoff en- 
ergy of the parton spectrum is reduced by about 10%, 
increasing the jet total cross section by about 50%. 



C. Physical interpretation — p-p 

The hard component from 200 GeV p-p collisions is 
well described by a parton spectrum folded with p-p frag- 
mentation functions measured at FNAL. The FF sup- 
pression at smaller pt /yt appears to be a real loss of part 
of the jet (the base), albeit a small fraction of the parton 
energy. Where is the missing jet base? 

For dijets in e-e collisions the color field is localized and 
continuous along the q-q axis. The fragment density on 
(from momentum component pz) along the color-dipole 
axis is uniform near the parton CM [46|. Correspond- 
ing FFs on y (from total momentum p) as in Fig. [3] (left 
panel) fall smoothly to zero at j/mm (pt ~ 0.05 GeV/c). 
In p-p collisions the "bases" of p-p FFs are missing rela- 
tive to e-e FFs, implying that the fragment distribution 
on yz acquires a "hole" at the parton CM. The hole in 
the dipole distribution suggests that the scattered-parton 
pair is not color connected, leading to these observations: 

1) Parton "scattering" in p-p collisions must proceed 
by exchange of a color singlet between nucleons (a variant 
of Pomeron exchange). All momentum transfers between 
nucleons are colorless, consistent with no color connec- 
tion between scattered partons. 

2) A scattered parton remains color connected to the 
parent nucleon. Part of the color field deviates from the 
parton-parton axis (follows the projectile nucleon). 

3) Fragmentation responds to the color connection, 
and some fragments are shifted away from the parton- 
parton axis (away from p-p mid-rapidity). FF bases are 
then missing from observed jets and hard component. 



The color-connection mechanism may explain suppres- 
sion of the bases of p-p FFs. If so, the larger the x of the 
scattered parton the larger should be the suppression or 
loss from the FF, as suggested by CDF FFs [HliS- 



D. Physical interpretation — Au-Au 

In more-central Au-Au collisions (above the sharp 
transition) the FD approaches a limiting case correspond- 
ing to "medium-modified" e-e FFs. Instead of suppres- 
sion at smaller yt there is enhancement. Fragmentation 
is modified, and jet bases are partially restored. The 
suppression mechanism proposed in the previous subsec- 
tion implies that color connection in more-central Au-Au 
collisions transitions from the N-N case toward the color 
dipole of e+-e" collisions. However, the splitting cascade 
is also modified within the Au-Au collision context. 

Minijet angular correlations reveal that restored jet 
bases are strongly elongated on rj relative to the larger-y^ 
part of the FF [3, [l3|- Triggered jet analysis refers to FF 
components as "jet" and "ridge" [i^. From this anal- 
ysis the 77 elongation can be explained by residual color 
connection to parent nucleons and resulting distortion of 
fragment distributions relative to the back-to-back dijet 
configuration in e^-e~ collisions. The larger the x of the 
struck parton the stronger should be the 77 elongation. 

The underlying parton spectrum remains the same ex- 
cept for the low-energy cutoff, which drops from 3 GeV 
to 2.7 GeV, increasing the jet cross section by 50%. Com- 
bined with a 3-fold increase in minijet multiplicity and 
doubling of dijet efficiency e the total fragment multiplic- 
ity per unit rj per N-N collision increases by a factor 9 in 
central Au-Au compared to p-p collisions. The spectrum 
cutoff may be reduced in more-central Au-Au collisions 
because the fragmentation process is altered: the hadron 
density of states becomes larger and the effective cutoff 
can then move down. There is no suggestion of parton 
loss by absorption or multiple scattering in a medium. 



E. Implications for conventional energy-loss models 

In a summary of high-pt physics at RHIC and LHC [i^ 
it was reported that "observation of the strong suppres- 
sion of high Pt hadrons" motivates the choice of 'hard 
probes" (energetic partons) to "characterize the decon- 
fined medium" at RHIC. Energetic partons "are ex- 
pected to lose energy through coUisional energy loss and 
medium- induced gluon radiation... known as jet quench- 
ing." Effects of a flowing medium may also be re- 
vealed jl^. But energy-loss schemes based on random 
multiple scattering of a leading parton are falsified by 
minijet angular-correlation data, especially reduction of 
the minijet azimuth width with increasing A- A centrality 
and survival of essentially all scattered partons above 3 
GeV as jet correlations in the hadronic final state. 
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Earlier high-pt results from RHIC such as "disappear- 
ance of the away-side jet" from dihadron correlations [H^l 
and "jet quenching" from spectra [5ll | were interpreted 
to imply absorption of most jets (including minijets) in 
central collisions within an "opaque core," observed jets 
being restricted to production in a surface layer. Those 
conclusions result in part from limitations of the mea- 
sures employed. From this analysis we find that all jets 
survive to particle detection, albeit some are modified. 

A measured negative shift of the entire FD on yt corre- 
sponds at most to a 25% reduction of the leading-parton 
energy (for central Au-Au). Because the parton spec- 
trum falls rapidly the shift does correspond to a 5-fold 
reduction of the FD at larger pt. But energy is not nec- 
essarily lost from integrated FFs, and the total hadron 
multiplicity is not suppressed, instead it is shifted to a 
different part of the pt spectrum and actually increases. 

Jet estimation by "leading particle" biases the jet 
structure (and imposes strong limitations on statistical 
power and accurate model tests), thus motivating full 
event-wise jet reconstruction. Regarding event-wise re- 
construction "...'jets' at RHIC are very complicated ob- 
jects which make them impossible to disentangle from the 
'background' " [emphasis added]. That rnay be true for 
event- wise reconstruction (except cf. [12,1131): but unbi- 
ased jet reconstruction is in effect accomplished on a com- 
binatoric basis via hard components of nuclear spectra 
(FDs) and minijet correlations. The observed structure 
can be unfolded to reveal the effective parton spectrum 
and FF ensemble as modified in nuclear collisions. 

"High-pt particle production in proton-proton colli- 
sions" is said to "provide the baseline 'vacuum' refer- 
ence to heavy-ion collisions to study the QCD medium 
properties" [4^|. But in this analysis we find that mis- 
taking the p-p hard component (or worse, the entire p-p 
spectrum) for an FD reference can distort "medium" ef- 
fects at smaller pt- Comparison of spectra to NLO FDs 
only above 2-3 GeV/c can also be misleading. The most 
significant modifications to fragmentation lie below that 
interval. Ironically, "soft probes" (minijets) reveal the 
most important details of QCD collision dynamics. 



F. Implications for hadrochemistry 

Strangeness and heavier flavors require exceptional en- 
ergy densities for exceptional production (relative to a 
statistical- model reference). Exceptional energy densities 
are accessible within minijets. Systematic strangeness 
trends, described as "strangeness suppression" in p-p 
collisions |54] and "strangeness enhancement" in more- 
central A-A collisions [55], suggest that minijet produc- 
tion may influence heavy-flavor (s, c, b) abundances in 
nuclear collisions. Because FF bases are suppressed in 
p-p collisions strangeness production is suppressed as 
well. In more-central A-A collisions the FF base is re- 
stored (even enhanced), and strangeness (and heavier fla- 
vors) are enhanced as well. In that scenario strangeness 



production would closely follow minijet production and 
"energy-loss" trends. A possible correspondence of mini- 
jet systematics with heavy-flavor abundances should be 
carefully considered before invoking the context of a ther- 
malized bulk medium, i.e., canonical suppression and en- 
hancement within an equilibrated statistical ensemble. 



G. Implications for saturation-scale arguments 

The saturation- scale model (SSM) (e.g. |4^]) is a 
limiting case of the two-component model of nuclear 
collisions in which the soft component disappears, all 
particle /pt/iJt production proceeds via the hard compo- 
nent, and all antecedent scattered partons are thermal- 
ized prior to hadronization. No predictions are made for 
hadronic correlations, and no justification is given for as- 
sumed disappearance of the soft component. 

According to the SSM hypothesis the parton spec- 
trum from A-A collisions is said to saturate at a par- 
ticular momentum scale po (saturation scale) depending 
on CM energy and A. Spectrum saturation is attributed 
to saturation of the initial parton density in the projec- 
tiles. Because of the parton spectrum structure almost 
all particle / Pt / Et production should then correspond to 
the saturation scale. The SSM parton spectrum cutoff for 
RHIC Au-Au colhsions is estimated to be ~ 1 GeV [ii ]. 

Semihard parton scatters from elementary hadronic 
collisions are isolated (<C one per collision), and p-p data 
imply that the density of hadronic final states constrains 
parton scattering near midrapidity. There is a direct QM 
coupling between parton scattering and hadronization in 
nuclear collisions: the latter constrains the former. 

In contrast, according to SSM any parton scatter that 
can happen kinematically will happen, independent of 
the density of hadronic final states, because there is no 
direct coupling between parton scattering and hadroniza- 
tion. Partons scatter to an intermediate QCD state (ex- 
panding QCD medium, possibly thermalized) from which 
hadrons later emerge in a collective process. The only 
constraint on parton scattering is the saturation limit of 
the initial parton flux density from the projectiles. 

Based on that reasoning the SSM parton spectrum is 
extended from 3 GeV down to 1 GeV, as illustrated in 
Fig. [7] (left panel) , implying a large increase in scattered 
partons compared to isolated p-p collisions (e.g., 85 mb 
dijet cross section in central Au-Au collisions vs 2.5 mb). 
The final-state ha,dron/ Et/pt production in A-A colli- 
sions is then attributed entirely to the hard component. 
The large soft-component production is incorrectly at- 
tributed to the SSM-invoked parton spectrum extension 
below 3 GeV to buttress the saturation-scale argument. 

The SSM hypothesis is contradicted by two-component 
analysis of correlations and spectra, particularly the cen- 
trality dependence (as in the present analysis). As 
demonstrated in Sec. lXI Bl the soft component (projectile 
nucleon fragmentation) dominates hadiOTi/pt/Et produc- 
tion even in central Au-Au collisions. The complemen- 
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tary hard-component production matches the ~ 3 GeV 
parton cutoff mferrcd from hard-component spectra. 

XIV. SUMMARY 

In this analysis fragment distributions measured as pt 
spectrum hard components are used to resolve and define 
both parton spectra and FF ensembles over large kine- 
matic intervals. Parton fragmentation and "energy loss" 
in nuclear collisions are studied. The initial conditions of 
nuclear collisions, particularly the scattered-parton en- 
ergy distribution relevant to pQCD and possibly hydro- 
dynamics, are estimated with reduced ambiguity. 

Parametrized p-p and e~^-e~ fragmentation functions 
are folded with a parton spectrum model to produce 
fragment distributions compared to hard components of 
hadron spectra. Comparison of a calculated p-p FD to 
a measured p-p hard component determines a reference 
parton spectrum for all nuclear collisions. The inferred 
spectrum agrees quantitatively with pQCD calculations. 

A theoretical model of FF medium modification is im- 
plemented by simple alteration of FF parametrizations. 
In-medium p-p and e~^-e~ FFs are folded with the refer- 
ence parton spectrum to produce modified FDs in turn 
compared to spectrum hard components for several cen- 
tralities of 200 GeV Au-Au collisions. Below a transition 
centrality hard components for peripheral Au-Au and p- 
p collisions are well-described by the in- vacuum p-p FD. 
Above the transition centrality hard-component data are 
well-described by in-medium e^-e~ FDs over the entire 
Pt spectrum [0.3,10] GeV/c. 

The implications of the analysis are as follows: 1) The 
underlying scattered-parton spectrum changes little from 
p-p to central Au-Au collisions; the exception is reduction 
of the 3 GeV cutoff energy by 10%. All partons survive 



to final-state manifestations, are not "thermalized." 2) 
Fragmentation functions for p-p collisions are substan- 
tially different from those for e^-e~ collisions — the low- 
momentum base is suppressed in the former, compris- 
ing 30-70% of the expected fragment number but only 
a few percent of the parton energy. 3) The missing FF 
base suggests that scattered partons from p-p collisions 
are not color connected to each other, that a scattered 
parton remains color connected to the parent projectile 
nucleon. 4) In more-central Au-Au collisions there is 
a sharp transition to in-medium e+-e~ FDs, suggesting 
that the color-connection topology changes substantially. 

5) Central Au-Au collisions appear to be nearly trans- 
parent to energetic partons and even to hadron fragments 
with Pt ^ 0.3 GeV/c which remain correlated with the 
parent parton, are not rescattered. 6) A measured large 
increase in Au-Au minijet angular correlations above the 
sharp transition on centrality corresponds to the large 
increase in mean dijet multiplicity for in-medium e"'"-e~ 
FFs compared to in-vacuum p-p FFs inferred from this 
analysis. 7) The most significant alteration of parton 
fragmentation in nuclear collisions occurs below pt = 2 
GeV/c. 8) A proposed extension of the scattered-parton 
Pt spectrum down to 1 GeV in central Au-Au collisions 
motivated by saturation-scale arguments is contradicted 
by spectrum and correlation data. 

Nuclear collisions have been described quantitatively 
in terms of perturbative QCD. Nonperturbative aspects 
already present in p-p collisions (fragmentation function 
modification and parton spectrum termination) evolve 
with heavy ion collision centrality in intuitively reason- 
able ways, leading to better understanding of QCD dy- 
namics over an extended space-time volume. 
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